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Abstract. We describe the propagation of supernova shocks within the surrounding 
medium, which may be due to mass-loss from the progenitor star. The structure and density 
profile of the ejected material and surrounding medium are considered. Shock wave interac- 
tion with clouds and wind-bubbles, and issues relevant to cosmic rays are briefly discussed. 
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1. Introduction 

The propagation of Supernova Remnants 
(SNRs) in the ambient medium has been stud- 
ied for over 4 decades. Early work was re- 
viewed by Woltjer (1972). Since then there 
have been several excellent reviews (Chevalier 
1977; Ostriker & McKee 1988; Chevalier & 
Fransson 1994; Bisnovatyi-Kogan & Silich 
1995). The high energy emission from SNRs 
was reviewed by Reynolds (2008). In this short 
presentation, our aim is to review the basic 
ideas of the propagation of SN shock waves 
into the ambient medium, and point out issues 
relevant to cosmic -ray acceleration. For further 
details we refer the reader to the excellent pub- 
lications cited above. 

2. Ejecta-Dominated Stage 

The explosion of massive stars (<: 8M ), or the 
thermonuclear deflagration and detonation of 
white dwarfs, leads to the production of highly 
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supersonic shocks waves that propagate into 
the ambient medium. In order to understand 
the evolution of these shock waves, one must 
know the density structure of the ejecta and 
the surrounding medium. The density struc- 
ture of the ejecta in core-collapse SNe can 
be approximated by a decreasing power law 
(Chevalier & Soker 1989; Matzner & McKee 
1999) with an index roughly between 9 and 
1 1 . Below a certain velocity the density is as- 
sumed roughly constant. The ambient medium 
in core-collapse SNe is much more compli- 
cated. Close to the star, the density can be ap- 
proximated by a wind, whose density varies 
as r~ 2 if the wind parameters are constant (al- 
though there is no reason for them to be con- 
stant, see Dwarkadas & Gruszko (2012)). The 
expansion of power-law ejecta decreasing as 
r~" into a power-law circumstellar medium 
whose density decreases as r _s can be de- 
scribed by a self-similar solution (Chevalier 
1982; Truelove & McKee 1999), where the ra- 
dius goes as t ( "~ 3) / ( "~' !) ~ f, where m < 1 
is called the expansion parameter. For repre- 
sentative values n — 9 and s — 2 indicative 
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of a wind, m = 0.86. Thus, the presence of 
a medium, however tenuous, restricts free ex- 
pansion of the shock wave, whose velocity de- 
creases as f 1 . Further from the star, the den- 
sity depends on the mass-loss from, and evo- 
lution of, the progenitor star, and may form a 
wind-blown bubble (Weaver et al. 1977, §5.2). 

Type la SNe are assumed to arise from 
low mass stars, which don't suffer from con- 
siderable mass-loss and don't modify their 
ambient medium significantly (Badenes et al. 
2007). A first assumption is to assume an 
unmodified constant density medium around 
la's. The ejecta profile is more complicated. 
By comparing density profiles obtained by hy- 
drodynamical simulations of Type la explo- 
sions, Dwarkadas & Chevalier (1998, hereafter 
DC98) showed that an exponential profile is 
a better approximation for Type la ejecta pro- 
files. The introduction of an exponential intro- 
duces one more variable, and therefore the so- 
lution is no longer self-similar, although it is 
still scalable, and can be expressed in terms of 
a single family of parameters. 

The expansion of the SN ejecta into the am- 
bient medium drives a shock into the medium. 
The material swept up by the shock is slower 
moving, and the ejecta need to be decelerated 
before they interact with this material. This 
leads to the formation of a reverse shock, that 
expands back into the ejecta in a Lagrangian 
sense. The two are separated by a contact dis- 
continuity, which divides the region of shocked 
ejecta from the region of shocked ambient 
medium (Fig. 1). The shocks will accelerate 
particles to high energies due to a process 
that is generally thought to be some variant of 
Diffusive Shock Acceleration. 

The self-similar, or scalable solutions, pro- 
vide the ratio of radii of the reverse shock, 
contact discontinuity (CD) and forward shock. 
These ratios are fixed for purely non-radiative 
hydrodynamical evolution, but may be modi- 
fied if significant energy (> 10%) is expended 
in accelerating particles to relativistic speeds 
(Caprioli 2011). Large modification of these 
ratios is often taken as evidence of cosmic -ray 
acceleration, which would tend to narrow the 
distance between the CD and forward shock 
for instance (Warren et al. 2005; Kosenko et al. 
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Fig. 1. The density (red) and pressure (purple) pro- 
files of the shocked interaction region between the 
forward and reverse shocks, for [Top] a Type la SN 
with an exponential ejecta density profile interact- 
ing with a constant density medium and [Bottom] A 
core-collapse SN with a power-law density profile 
interacting with a wind with constant parameters. 
Time in years is given at the top. Note the variation 
in the profile between the two cases. An expanding 
grid is used to allow most of the computational do- 
main to be occupied by the shocked region. 

2011). However this presumes that the ejecta 
profiles, and especially the smooth, homoge- 
neous ambient medium profiles, are exact rep- 
resentations. For Type la's, DC98 showed that 
while the exponential is a good and conve- 
nient approximation to the ejecta profiles, de- 
viations from this are important, especially for 
low-energy explosions. One would also ex- 
pect inhomogeneities in the ambient medium 
over large scales, especially for core-collapse 
SNe, whose effect may remain imprinted in the 
profile over several doubling times, and thus 
modify the ratio. Other factors, including the 
growth of R-T instabilities (Orlando, these pro- 



Dwarkadas: SNe in CSM and ISM 



783 



ceedings), could alter the distance between CD 
and forward shock (Kosenko et al. 201 1). 

Since the reverse shock is generally ex- 
panding into a higher density than the forward 
shock, it's velocity is lower. If particles are ac- 
celerated by the reverse shock, as suggested 
for Cas A (Helder & Vink 2008), they would 
be expected to have a lower maximum en- 
ergy than those at the forward shock. If core- 
collapse SNe expand in winds whose density 
decreases as r~ 2 , whereas Type la SNe expand 
into a constant density medium, then the ra- 
tio of the density ahead of the reverse shock 
to that ahead of the forward shock falls more 
rapidly in Type la's then in core-collapse SNe, 
suggesting that the reverse shock's effects are 
felt longer in core-collapse SNe than in la's. 

During the early stages of the evolution, 
both the forward and reverse shocks expand 
outwards in radius. The decelerating contact 
discontinuity is Rayleigh-Taylor (R-T) unsta- 
ble, and R-T fingers will mix shocked ejecta 
material with the shocked ambient medium 
(Chevalier et al. 1992; Dwarkadas 2000). At 
a later time, depending on the ejecta profile 
(Dwarkadas & Chevalier 1998), the ejecta den- 
sity becomes low enough that the reverse shock 
will "turn over" and actually begin to move in- 
wards in radius, towards the explosion center. 
It is only after the reverse shock reaches the 
center that the SNR can be thought of as enter- 
ing the Sedov-Taylor phase. 

3. Sedov-Taylor Stage 

In this adiabatic stage the remnant expansion 
depends mainly on two parameters, the energy 
of the explosion E and the ambient medium 
density p. The shock front travels as R s h °c 
(E/p) 1 ^ 5 1 1 ! 5 when expanding into a constant 
density medium, or R s j, cc t 2 ^ in the less likely 
case of the stage being reached while the SN 
is expanding in a wind medium with p oc r 2 . 
This expression was first derived for a point ex- 
plosion in a uniform medium (Taylor 1950). 

In order for the remnant to reach the Sedov 
stage, several things must happen (1) The re- 
verse shock must reach the center and effec- 
tively cease to exist, since the Sedov-Taylor so- 
lution is a single-shock solution, as opposed to 



the double-shocked solutions earlier (Fig. 1). 
In multi-dimensional simulations, the reverse 
shock generally does not remain spherical as it 
moves back, and eventually tends to dissipate. 
(2) The expansion parameter m discussed ear- 
lier must decrease from its larger value till it 
reaches the Sedov value of 0.4 (0.67) in a con- 
stant density medium (wind) (3) The density 
profiles must change from those of the ejecta- 
dominated stage (Fig. 1) to those of the Sedov 
stage. Although traditionally the Sedov stage 
is assumed to be reached when the swept-up 
mass is just larger than the ejecta mass, this 
does not hold in practice, as first shown almost 
4 decades ago (Gull 1973). DC98 showed that 
the swept-up mass must be about 20-30 times 
the ejecta mass (depending on the ejecta pro- 
file) when the reverse shock reaches the center. 

Thus the idealistic representation of the 
Sedov stage is reached later in time, and SN ra- 
dius, than when the swept-up mass equals the 
ejecta mass. The simple reason is that a SNR 
is not a point explosion, and the Sedov stage 
arises from the ejecta-dominated stage, when 
the remnant is already light-years across in ra- 
dius. Major galactic remnants, such as Tycho, 
Cas A, Kepler, and even the 1000-year old 
SN 1006 are currently in the ejecta-dominated 
stage, with a double-shocked structure. 

Some calculations suggest that the trans- 
formation of SNR energy to cosmic -rays is 
most efficient at the beginning of the Sedov 
stage (Berezhko & Volk 1997; Ptuskin & 
Zirakashvili 2005; Ptuskin et al. 2010). This 
has prompted many investigations into the ac- 
celeration of particles in the Sedov stage (Kang 
2006, 2010; Castro et al. 201 1). The transition 
into the Sedov stage however may occupy a 
significant amount of the lifetime of this stage, 
and must be taken into account in computa- 
tions of the accelerated particle spectrum. In 
some cases, the Sedov stage may be shortened, 
or even completely eliminated (§5.2). 

4. Radiative Stage: 

When the shock slows down sufficiently that 
the cooling time of the material behind the 
shock becomes smaller than the flow time, the 
SNR enters the radiative phase. This phase 
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is not generally of much interest in cosmic- 
ray physics because the slow shock is inef- 
ficient at accelerating particles to relativistic 
energies (Berezhko & Volk 1997). The evo- 
lution can be approximated by a momentum- 
driven flow, where R s i t oc t 1/4 . Ostriker & 
McKee (1988) prefer a pressure-driven snow- 
plow model (R s i, oc f 2 ^ 7 ), while Cioffi et al. 
(1988) suggest an offset solution R s h oc t 03 . 
In any case the SNR has decelerated consider- 
ably from the initial stages. Eventually, it will 
merge with the ISM, finally mixing all the con- 
stituents of the explosion back into the ISM. 

A radiative shock though can form at any 
time depending on the density of the ambient 
medium, if the cooling time t coo i oc lcT s /(nA) 
becomes smaller than the flow time. Here T s oc 
v] is the postshock temperature, n the den- 
sity of the ambient medium, and A the cool- 
ing function. This can happen for young SNe 
expanding into a high density medium, as has 
been postulated for SN 1993 J (Nymark et al. 
2009; Chandra et al. 2009). SNRs impacting 
with clouds or driven into clumps may also 
form radiative shocks. For some SNe, mod- 
elling of the X-ray and optical emission indi- 
cates a radiative shock impacting clumps while 
an adiabatic shock expands into the interclump 
medium (Chugai 1993; Chugai & Danziger 
1994; Chugai et al. 1995). 

5. Some special Case Studies 

5.1. SNR-Molecular Cloud Interaction 

This is an interesting case for cosmic ray physi- 
cists, because many SNRs discovered in y- 
rays are found to be interacting with molec- 
ular clouds, suggesting a hadronic origin for 
the emission. There are various possibilities, 
depending on whether the SNR shock is ac- 
tually interacting with the cloud or not. (1) 
Shock-cloud interaction has not occurred. The 
escaped particles from the SNR have inter- 
acted with the dense cloud material, leading to 
broadband non-thermal emission (Gabici et al. 
2009). The hydrodynamical evolution of the 
SNR is not yet affected by the presence of 
the cloud. (2) The SN shock has impacted the 
cloud. The impact of SNR shocks with clouds 



was discussed by Klein et al. (1994), and re- 
sulting y-ray emission by Uchiyama et al. 
(2010). Fig 2 shows results from a simula- 
tion of a SNR interacting with a cloud of den- 
sity pci. The impact drives a shock into the 
cloud, and a reflected shock back into the in- 
tercloud medium with density p !C . The velocity 
of the cloud shock v c ; oc v (c V(p, e /Pd) can be 
much smaller than the inter-cloud shock veloc- 
ity. In case of a dense cloud, the cloud shock 
can become radiative if it cools faster than the 
flow time through the cloud. The impact of 
the shock on the cloud, and the shear flow on 
the front and sides, leads to the formation of 
Richtmeyer-Meshkov and Kelvin-Helmholtz 
instabilities, which will eventually lead to the 
destruction of the cloud. On a smaller scale, 
these types of interactions are seen in SNRs in- 
teracting with a clumpy medium. 

5.2. SNRs in Wind-blown Bubbles 

Many massive stars will form large, and some 
not-so-large, wind bubbles around the star. 
Wind bubbles are regularly seen around main- 
sequence and Wolf-Rayet (W-R) stars (Chu 
2008). The basic structure (Fig 3) comprises of 
an inner and outer shock separated by a contact 
discontinuity (CD). The shocked wind between 
the inner shock and CD usually has a low den- 
sity and high pressure, and thus high tempera- 
ture. The outer shock is usually radiative, and 
the material between the outer shock and CD 
forms a thin, dense shell. Ionization from the 
star may result in the formation of an ionized 
region around the wind-blown shell. 

The post-main-sequence phases occupy a 
small fraction of the total stellar lifetime, but 
may result in substantial mass-loss and modifi- 
cation of the ambient medium (van Marie et al. 
2005; Freyer et al. 2006; Dwarkadas 2007b; 
Toala & Arthur 2011). Stellar evolution the- 
ory suggests that single stars < 30 M© will end 
their lives as RSGs to form Type IIP SNe. Stars 
with higher initial mass will become W-R stars, 
with perhaps an intermediate Luminous Blue 
Variable phase, and explode as Type Ib/c SNe. 
RSGs have slow dense winds. W-R stars have 
fast winds with a lower mass-loss rate (Puis, 
these proceedings), resulting in a lower den- 
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Fig. 2. Snapshots from a simulation of the interaction of a SNR with a dense cloud, carried out using the 
VH-1 code in spherical symmetry, (a) A SNR with a spherical double-shocked structure approaches the 
cloud. Time in years is at top of picture, (b) The impact drives a shock into the cloud, and a reflected 
shock back into the inter-cloud medium. A complicated hydrodynamic structure is seen, (c) The rest of 
the intercloud shock has swept past, while the reflected shock can be clearly seen. Meanwhile, the shock 
crossing the cloud results in the formation of instabilities that will ultimately destroy the cloud. 
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Fig. 3. The basic structure of a wind-blown bub- 
ble around a main-sequence 40 M Q star, carried 
out using the ionization-gasdynamics code AVATAR 
(Dwarkadas & Rosenberg 2011). Density in red 
(scale on left), ionization fraction in blue (scale on 
right). Radius in parsecs on X-axis. From left, we 
see the freely expanding wind, the wind termination 
shock, a shocked wind region, an ionized HII region, 
and the dense shell. The structure upto the HII re- 
gion is fully ionized. The ionization front is trapped 
within the shell, and the ionization fraction can be 
seen to drop very sharply within the shell. 



sity in the ambient medium into which the SN 
will expand. A binary companion can alter the 
stellar evolution and structure of the ambient 
medium, as postulated for SN 1987 A (Morris 
& Podsiadlowski 2007). 

The evolution of SN shock waves 
within these bubbles has been studied by 



Chevalier & Liang (1989); Tenorio-Tagle et al. 
(1990, 1991); Dwarkadas (2005, 2007a,b); 
Dwarkadas et al. (2010). The collision of a 
SN shock with a shock or CD results in a 
transmitted shock expanding outwards, and a 
reflected shock expanding back into the ejecta. 
Thus expansion within a bubble results in 
several shocks bouncing around in the remnant 
cavity. Many of these shocks can accelerate 
particles to high energies, at least for short 
periods of time. Accelerated particles that 
escape from the SNR and interact with this 
shell may give rise to y-ray emission (Ellison 
& Bykov 2011). When the SN shock collides 
with the dense shell, it can result in significant 
deceleration of the shock, and an increase in 
the optical and X-ray emission, as seen in 
SN 1987A (Park et al. 2006) and SN 1996cr 
(Bauer et al. 2008). The shock will traverse 
the shell with a low velocity (see 5.1) until it 
emerges. In extreme cases, the shell density 
may be so high that the shock will be trapped 
in the shell, and the kinetic energy may be 
radiated away. The SN will then go directly to 
the radiative phase, avoiding the Sedov phase. 

Acknowledgements. I am very grateful to the orga- 
nizers, especially A. Marcowith, for organizing a 
splendid conference, for inviting me to present this 
review, and for their gracious hospitality. It is a plea- 
sure to acknowledge the many people from whom 
I have learned about SNe and SNRs, especially R. 
Chevalier. My research is supported by grants from 
Chandra, and Fermi grant NNX10AO44G. 



786 



Dwarkadas: SNe in CSM and ISM 



References 

Badenes, C, Hughes, J. P., Bravo, E., & 

Langer, N. 2007, ApJ, 662, 472 
Bauer, F. E., Dwarkadas, V. V., Brandt, W. N., 

etal. 2008, ApJ, 688, 1210 
Berezhko, E. G. & Volk, H. J. 1997, 

Astroparticle Physics, 7, 183 
Bisnovatyi-Kogan, G. S. & Silich, S. A. 1995, 

Reviews of Modern Physics, 67, 661 
Caprioli, D. 201 1, JCaP, 5, 26 
Castro, D., Slane, P., Patnaude, D. J., & 

Ellison, D. C. 201 1, ApJ, 734, 85 
Chandra, P., Dwarkadas, V. V., Ray, A., 

Immler, S., & Pooley, D. 2009, ApJ, 699, 

388 

Chevalier, R. A. 1977, ARA&A, 15, 175 
Chevalier, R. A. 1982, ApJ, 258, 790 
Chevalier, R. A., Blondin, J. M., & Emmering, 

R.T. 1992, ApJ, 392, 118 
Chevalier, R. A. & Fransson, C. 1994, ApJ, 

420, 268 

Chevalier, R. A. & Liang, E. P. 1989, ApJ, 344, 
332 

Chevalier, R. A. & Soker, N. 1989, ApJ, 341, 
867 

Chu, Y. 2008, in IAU Symposium, Vol. 
250, IAU Symposium, ed. F. Bresolin, 
P. A. Crowther, & J. Puis, 341-354 
Chugai, N. N. 1993, ApJ, 414, L101 
Chugai, N. N. & Danziger, I. J. 1994, MNRAS, 
268, 173 

Chugai, N. N, Danziger, I. J., & della Valle, 

M. 1995, MNRAS, 276, 530 
Cioffi, D. F, McKee, C. F, & Bertschinger, E. 

1988, ApJ, 334, 252 
Dwarkadas, V. & Rosenberg, D. 2011, in 

BAAS, Vol. 43, American Astronomical 

Society Meeting Abstracts 217, 251.13 — h 
Dwarkadas, V. V. 2000, ApJ, 541, 418 
Dwarkadas, V. V. 2005, ApJ, 630, 892 
Dwarkadas, V. V. 2007a, Ap&SS, 307, 153 
Dwarkadas, V. V. 2007b, ApJ, 667, 226 
Dwarkadas, V. V. & Chevalier, R. A. 1998, 

ApJ, 497, 807 
Dwarkadas, V. V., Dewey, D., & Bauer, F. 

2010, MNRAS, 407, 812 
Dwarkadas, V. V. & Gruszko, J. 2012, 

MNRAS, 419, 1515 
Ellison, D. C. & Bykov, A. M. 201 1, ApJ, 731, 



87 

Freyer, T., Hensler, G, & Yorke, H. W. 2006, 

ApJ, 638, 262 
Gabici, S., Aharonian, F. A., & Casanova, S. 

2009, MNRAS, 396, 1629 
Gull, S. F. 1973, MNRAS, 161, 47 
Helder, E. A. & Vink, J. 2008, ApJ, 686, 1094 
Kang, H. 2006, Journal of Korean 

Astronomical Society, 39, 95 
Kang, H. 2010, Journal of Korean 

Astronomical Society, 43, 25 
Klein, R. I., McKee, C. F, & Colella, P. 1994, 

ApJ, 420,213 
Kosenko, D., Blinnikov, S. I., & Vink, J. 201 1, 

ArXiv e-prints 
Matzner, C. D. & McKee, C. F. 1999, ApJ, 510, 

379 

Morris, T. & Podsiadlowski, P. 2007, Science, 
315, 1103 

Nymark, T. K, Chandra, P., & Fransson, C. 

2009, A&A, 494, 179 
Ostriker, J. P. & McKee, C. F. 1988, Reviews 

of Modern Physics, 60, 1 
Park, S., Zhekov, S. A., Burrows, D. N., et al. 

2006, ApJ, 646, 1001 
Ptuskin, V, Zirakashvili, V, & Seo, E.-S. 2010, 

ApJ, 718, 31 
Ptuskin, V. S. & Zirakashvili, V. N. 2005, 

A&A, 429, 755 
Reynolds, S. P. 2008, ARA&A, 46, 89 
Taylor, G. 1950, Royal Society of London 

Proceedings Series A, 201, 159 
Tenorio-Tagle, G, Bodenheimer, P., Franco, J., 

& Rozyczka, M. 1990, MNRAS, 244, 563 
Tenorio-Tagle, G, Rozyczka, M., Franco, J., & 

Bodenheimer, P. 1991, MNRAS, 251, 318 
Toala, J. A. & Arthur, S. J. 2011, ArXiv e- 

prints 

Truelove, J. K. & McKee, C. F. 1999, ApJS, 
120, 299 

Uchiyama, Y, Blandford, R. D., Funk, S., 
Tajima, H., & Tanaka, T. 2010, ApJ, 723, 
L122 

van Marie, A. J., Langer, N., & Garcfa-Segura, 

G. 2005, A&A, 444, 837 
Warren, J. S., Hughes, J. P., Badenes, C, et al. 

2005, ApJ, 634, 376 
Weaver, R., McCray, R., Castor, J., Shapiro, P., 

& Moore, R. 1977, ApJ, 218, 377 
Woltjer, L. 1972, ARA&A, 10, 129 



